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Abstract. We show that the Liquid Drop Model is best suited to describe the masses of 
prolate deformed nuclei than of spherical nuclei. To this end three Liquid Drop Mass formulas are 
employed to describe nuclear masses of eight sets of nuclei with similar quadrupole deformations. 
It is shown that they are able to fit the measured masses of prolate deformed nuclei with an 
RMS smaller than 750 keV, while for the spherical nuclei the RMS is, in the three cases, larger 
than 2000 keV. The RMS of the best fit of the masses of semi-magic nuclei is also larger than 
2000 keV. The parameters of the three models are studied, showing that the surface symmetry 
term is the one which varies the most from one group of nuclei to another. In one model, isospin 
dependent terms are also found to exhibit strong changes. The inclusion of shell effects allows 
for better fits, which continue to be better in the prolate deformed nuclei region. 



1. Introduction 

The description of nuclear masses in terms of the Liquid Drop Model paved the way to the basic 
understanding of nuclear properties, like the saturation of the nuclear force, the existence of 
pairing and shell effects, and the description of fission and fusion processes pQ. The Q-values of 
different nuclear reactions, obtained from mass differences, must be accurately known to allow 
the description of the astrophysical origin of the elements [2] . Accurate theoretical predictions 
of nuclear masses remain a challenge [3], sharing the difficulties with other quantum many-body 
calculations, and complicated by the absence of a full theory of the nuclear interaction. 

Decades of work have produced microscopic and macroscopic mass formulas [I]. At present, 
the most successful approaches seem to be the Finite Range Droplet Model (FRDM) [5], the 
Skyrme and Gogny Hartee Fock Bogolyubov (HFB) [61 [7], and the Duflo-Zuker (DZ) mass 
formula [H I10j . They allow for the calculation of masses, charge radii, deformations, and 
in some cases also fission barriers. They all contain a macroscopic sector which resembles the 
Liquid Drop Mass (LDM) formula, and include deformation effects. HFB calculations are now 
able to fit known nuclear masses with deviations competitive with the '95 FRDM calculations, 
which are also being improved, while the most precise and robust nuclear mass predictions are 
given by the DZ model [H 111]. 

The Liquid Drop Mass (LDM) formula captures the macroscopic features of the mass 
dependence on the number of neutrons N, of protons Z, and on its mass numbers A = N + Z. 



It includes volume and surface terms, the Coulomb interaction between protons, Wigner and 
symmetry terms, linear and quadratic in the neutron excess N — Z, and a pairing term. It is 
generally assumed that the liquid-drop energy of a spherical nucleus is described by a Bethe- 
Weizsacker mass formula |12| . being a common practice to describe nuclear masses and radii of 
spherical closed-shell nuclei in terms of a mean field and add deformation and other shell effects 
as corrections [T3] , 

It is the purpose of this contribution to show that, when nuclei with measured masses 
are grouped according with their quadrupole deformations, three Liquid Drop Mass formulas 
consistently allow for a fit with an RMS smaller than 750 keV for the set of most prolate deformed 
nuclei, while they are unable to fit the masses of spherical nuclei with an RMS smaller than 
2000 keV. Semi-magic nuclei are fit with a similarly large RMS. The parameters of the three 
models are studied, showing that the surface symmetry term is the one which varies the most 
from one group of nuclei to another. In one model, isospin dependent terms are also found to 
exhibit strong changes. 

Shell effects refer to the differences between the experimental binding energies |14| and the 
LDM predictions. In this work they are included employing linear, quadratic [15^ I16j . 3- and 
4-body terms [13], functions of the number of valence nucleons, following the ideas of the Duflo- 
Zuker model [9j [10]. The inclusion of shell effects allows for better fits, while the smallest rms 
are found again in the prolate deformed nuclei region. 



2. The three Liquid Drop Mass formulas 

Three Liquid Drop Mass formulas will be employed to analyze their ability to fit nuclear masses. 
One of them treats consistently volume and surface effects [13] . The second one incorporates 
explicitly isospin effects [12] while the third one includes quadratic isospin effects, the diffuseness 
correction to the Coulomb energy, the charge exchange correction term and the curvature energy 

In Ref. [T3] an improved version of the liquid-drop mass formula with modified symmetry 
and Coulomb terms is built, following a consistent treatment of nuclear bulk and surface effects. 
The negative nuclear interaction energy is given by: 

a 42/3 a 4T(T+1) Z(Z-l) A 

E LDM1 = -a v A + a s Ai + S v - ^_ 1/3) + a c ^ _ - a p ^, (1) 

where the pairing interaction is given by A = 2, 1, and for even-even, odd-mass and odd-odd 
nuclei, respectively. A modification A in the Coulomb term is included 

N-Z 5tt 2 d 2 

6Z(l + y- 1 A 1 /3) 6 rgA 2 / 3 ' 1 ' 

Its presence in the denominator of the Coulomb term suggests that it can be viewed as a 
correction to the radius of the nucleus. The symmetry term employs 4T(T + 1) instead of 
(N — Z) 2 to account for the Wigner energy. The Coulomb interaction is proportional to Z{Z — 1) 
to avoid the Coulomb interaction of a proton with itself. 

The second LDM formula we are going to analyze was introduced in Ref. [12] by considering 
isospin effects. In this case the liquid-drop energy of a spherical nucleus is described by a 
modified Bethe-Weizsacker mass formula 

E L DM2 = ~a v A + a s A 2 / 3 + a sym I 2 A + a c y /3 ' (1 - Z 2 ' 3 ) - a pair ^, (3) 



with isospin asymmetry / = (N — Z)/A. The pairing term is taken form [18J 
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based on the conventional surface-symmetry term of liquid-drop model, with a small correction 
term for description of isospin dependence of a sym . The symmetry energy coefficient d S y m 
increases with decreasing isospin asymmetry. This I correction term approximately describes 
the Wigner effect for heavy nuclei. 

Different mass formulae derived from the liquid drop model and the pairing and shell energies 
of the Thomas- Fermi model have been studied and compared in Ref [T7]. We selected for this 
study to include the diffuseness correction to the Coulomb energy, the charge exchange correction 
term and the curvature energy. In Ref. [T7] it is reported that the Coulomb diffuseness correction 
Z 2 /A term and the charge exchange correction Z 4 '/ 3 /A 1 / 3 term play the main role to improve 
the accuracy of the mass formula. The Wigner term and the curvature energy can also be used 
separately but their coefficients are very unstable. 

Their LDM formula is 

E LDM3 {N,Z)= -a v {l-k v I 2 )A + a s {l-k s I 2 )A^ + a k {l-k k I 2 )A^ + l^^ 
—fpZ 2 '/A — a CfiXC Z^I 3 1 A 1 / 3 + E pa i r 

The volume energy corresponding to the saturated exchange force and infinite nuclear matter 
is given by the first term. I 2 A is the asymmetry energy of the Bethe-Weizsacker mass formula. 
The second term is the surface energy. Its origin is the deficit of binding energy of the nucleons 
at the nuclear surface and corresponds to semi-infinite nuclear matter. The following term is 
the curvature energy. It results from non-uniform properties which correct the surface energy 
and depends on the mean local curvature. The decrease of binding energy due to the Coulomb 
repulsion is given by the fourth term, which has and adjustable charge radii r^A 1 ^ . The Z 2 /A 
term is the diffuseness correction to the sharp radius Coulomb energy, also called also the proton 
form-factor. The Z^/A 1 / 3 term is the charge exchange correction term. The pairing energies 
of the Thomas-Fermi model |19| were employed. 

3. The fits 

The coefficients of the three LDM were selected to minimize the root mean square deviation 
(RMS) when the predicted binding energies BE t h(N Z) are compared with the experimental 
ones BE exp (N, Z), reported in AME03 |14| . modified so as to include more realistically the 
electron binding energies as explained in Appendix A of Lunney, Pearson and Thibault [4J. 

RMS [ E lE^d^l Z ) ~ BE ^ N > Z )] 2 ] V2 , (7) 

N nuc i is the number of nuclei in each group, listed in the fourth row of Table [3l The minimization 
procedure uses the routine Minuit [20] . 

The fits were performed employing the masses of nine groups of nuclei: 



group 


all 


1 


2 


3 


4 


5 


6 


7 


semi-magic 


e2 min 


-0.65 


-0.65 


-0.11 


0.00 


0.04 


0.12 


0.18 


0.23 




e2 max 


0.65 


-0.11 


0.00 


0.04 


0.12 


0.18 


0.23 


0.65 






2149 


258 


252 


332 


272 


307 


364 


364 


185 



Table 1. The nine groups on nuclei employed in the present study, their range of quadrupole 
deformation, and their number of nuclei. 

• all nuclei whose measured masses are reported in AME03 [H], which have N, Z > 8, 

• seven groups of nuclei whose quadrupole deformations, taken form the FRDM [21], lie in 
the ranges listed in the second and third row of Table [3j 

• the group of all semi-magic nuclei, having Z= 8, 20, 28, 50, 82 or N= 8, 20, 28, 50, 82 or 
126. 

Notice that group 1 contains most of the oblate nuclei, that the more spherical nuclei belong to 
group 3, and that the more prolate deformed nuclei are included in groups 6 and 7. 

For each LDM equation, nine fits were performed, one for each group of nuclei. In this 
way, nine sets of parameters were obtained, which minimize the RMS of each group of nuclei. 
Employing these nine sets of parameters, the RMS were estimated for all groups, whose RMS 
are shown in the next subsections. 

3.1. Analysis of the LDM1 formula 

Here we show the results obtained using Eq. dl|) for the nine regions. In our calculations we 
select nine sets of fixed values for the parameters a v , a s , o c , a p , S v and y. The results are 
exhibited in Tables [2] and [3j 

Table 2. RMS (in keV) for the nine groups (columns), for the nine sets of parameters (rows), 
employing Eq. ([TJ 





all 


1 


2 


3 


4 


5 


6 


7 


semi-magic 


set a u 


2387 


1842 


3814 


3080 


2037 


1495 


1815 


2060 


3888 


seti 


3300 


1313 


3113 


3420 


3277 


3173 


5194 


1476 


4016 


set 2 


4053 


2917 


1676 


2433 


3088 


4273 


6243 


4701 


2570 


set 3 


3404 


3313 


2298 


2063 


2556 


3408 


4016 


4721 


2600 


set4 


2567 


2431 


3867 


3096 


1746 


1741 


1975 


2626 


3901 


sets 


2630 


2097 


4848 


3696 


2072 


1053 


1360 


1735 


4521 


set6 


2828 


1942 


5205 


4262 


2578 


1469 


870 


1296 


5078 


sety 


3169 


1655 


4077 


4220 


3543 


2753 


3643 


746 


5016 


set sem j 


4514 


3816 


1930 


2665 


3316 


4777 


6729 


5306 


2113 



3.2. Analysis of the LDM2 formula 

Here we show the results obtained using Eq. @ for the nine regions. In our calculations we 
select nine sets of fixed values for the parameters a v , o s , o c , (ipain Csym &nd k. The results are 
exhibited in Tables 0] and EJ 



Table 3. Sets of parameters which minimize the RMS for the nine groups of nuclei, employing 
Eq. (JTJ). In the last three rows the average value of each parameter, their dispersion and 
percentage variation %var = 100 ^era.^ are shown. 





a v 


a s 


a c 


a p 


S v 


y 


set a // 


15.822 


18.491 


0.703 


6.048 


30.701 


2.616 


seti 


15.515 


17.577 


0.686 


6.092 


25.997 


1.583 


set 2 


15.883 


18.622 


0.712 


5.733 


28.711 


2.284 


set 3 


15.978 


18.884 


0.716 


5.895 


32.091 


3.059 


set4 


15.879 


18.564 


0.708 


4.103 


31.459 


2.666 


sets 


15.832 


18.435 


0.707 


5.266 


30.575 


2.482 


set6 


15.609 


17.835 


0.689 


5.334 


29.762 


2.388 


set7 


15.553 


17.782 


0.684 


5.126 


27.402 


1.865 


Set^gmj 


15.775 


18.012 


0.712 


5.500 


26.407 


1.489 


average 


15.761 


18.245 


0.702 


5.455 


29.234 


2.270 


a 


0.153 


0.426 


0.012 


0.579 


2.093 


0.495 


% var 


1.0 


2.3 


1.7 


10.6 


7.1 


21.8 



Table 4. RMS (in keV) for the nine groups (columns), for the nine sets of parameters (rows), 
employing Eq. © 





all 


1 


2 


3 


4 


5 


6 


7 


semi-magic 


set a // 


2374 


1788 


3733 


3081 


2044 


1502 


1855 


2058 


3916 


seti 


3208 


1254 


3114 


3400 


3243 


3077 


4947 


1435 


4060 


set 2 


4092 


2965 


1675 


2380 


3027 


4318 


6393 


4697 


2380 


set 3 


3475 


3352 


2173 


2069 


2596 


3518 


4183 


4816 


2551 


set4 


2561 


2408 


3761 


3052 


1762 


1771 


2036 


2688 


3897 


sets 


2616 


2073 


4823 


3702 


2083 


1021 


1311 


1704 


4548 


set6 


2795 


1865 


5130 


4234 


2619 


1435 


838 


1242 


5090 


set7 


3074 


1589 


4216 


4238 


3541 


2590 


3138 


656 


5093 


set sem j 


4677 


3915 


1959 


2722 


3427 


4953 


7036 


5466 


2056 



3.3. Analysis of the LDM3 formula 

Here we show the results obtained using Eq. (0) for the nine regions. In our calculations we 
select nine sets of fixed values for the parameters a v , o s , ro, ap a i r , k v , k s , fp, ac,exa 

a k andA; fc . 

The results are exhibited in Tables [6] and [7J 
3.4- LDM and deformation 

Tables [21 [H [6] display the RMS obtained with the three LDM formulas. Each row refers to one 
fixed set of parameters, each column to one group of nuclei. 

For the LDM1, Table El the global RMS is 2.39 MeV for the 2149 nuclei. In some groups, 
containing around two or three hundred nuclei, the RMS obtained with this set of parameters 
is smaller than the global one. The largest RMS are found in groups 2 and 3, those containing 
spherical nuclei. Along each column, corresponding to one group of nuclei, the smallest RMS 



Table 5. Sets of parameters which minimize the RMS for the nine groups of nuclei, employing 
Eq. ([3]). In the last three rows the average value of each parameter, their dispersion and 
percentage variation %var = 100 ^era.^ are snown - 





a v 


a s 


a c 


O'pair 


C-sym 


K 


set a // 


15.711 


18.920 


0.720 


6.989 


30.045 


1.587 


seti 


15.432 


18.119 


0.700 


6.793 


26.488 


1.260 


set 2 


15.695 


18.802 


0.722 


7.443 


27.979 


1.441 


set 3 


15.903 


19.395 


0.737 


6.717 


30.582 


1.642 


set4 


15.729 


18.877 


0.721 


4.685 


30.858 


1.635 


sets 


15.685 


18.756 


0.721 


6.210 


29.969 


1.548 


set6 


15.493 


18.254 


0.703 


6.276 


29.490 


1.553 


set7 


15.500 


18.398 


0.701 


6.193 


27.715 


1.361 


Set^gmj 


15.627 


18.329 


0.724 


5.327 


26.043 


1.099 


average 


15.642 


18.650 


0.717 


6.293 


28.797 


1.458 


a 


0.138 


0.383 


0.012 


0.801 


1.687 


0.175 


% var 


0.9 


2.0 


1.7 


12.7 


5.8 


12.0 



Table 6. RMS (in keV) for the nine groups (columns), for the nine sets of parameters (rows), 
employing Eq. ^ 





all 


1 


2 


3 


4 


5 


6 


7 


semi-magic 


set a // 


2422 


1884 


3773 


3186 


2031 


1519 


1861 


2125 


3937 


seti 


3576 


1183 


3067 


3785 


3461 


3456 


5854 


1283 


4424 


set 2 


4211 


3032 


1597 


2648 


2928 


4409 


6596 


4871 


2333 


set 3 


3563 


3398 


2197 


2151 


2618 


3628 


4455 


4808 


2499 


set4 


3049 


3035 


3849 


4090 


1517 


1764 


2560 


3456 


3952 


sets 


2776 


2406 


4838 


4004 


1927 


986 


1587 


2110 


4464 


set6 


2895 


2083 


5073 


4459 


2701 


1471 


819 


1579 


5103 


set7 


3200 


1484 


4128 


4317 


3506 


2701 


3790 


629 


5135 


set sem j 


4707 


3776 


1781 


2720 


3337 


4962 


7321 


5392 


1967 



always corresponds the set of parameters obtained fitting in this group, as expected. These RMS 
values are displayed in bold numbers. It is remarkable that the smallest RMS values are found 
in the two groups having the more oblate deformed nuclei, in groups 6 and 7, with 0.87 and 
0.75 MeV, respectively. On the other hand, the best fit of the spherical nuclei in group 3 has an 
RMS larger than 2.0 MeV, and those nuclei with very small quadrupole deformation, belonging 
to groups 2 and 4, have RMS larger than 1.6 MeV. Consistently, the semi-magic group of nuclei, 
which have small quadrupole deformations, have an RMS of 2.1 MeV. 

By analyzing the RMS for each column, it is possible to notice that the smaller RMS are 
always found around the smallest one. It supports the idea that the division in the seven groups 
with different deformations makes sense, because the parameters obtained fitting nuclei with 
a close deformation produce also a small RMS. Notice, for example, that among the RMS of 
groups 1 and 7, the more oblate and prolate deformed, respectively, the largest RMS are found 
with the sets of parameters 2 and 3, i.e. those fitted for spherical nuclei. On the opposite 



Table 7. Sets of parameters which minimize the RMS for the nine groups of nuclei, employing 
Eq. ([6]). In the last three rows the average value of each parameter, their dispersion and 
percentage variation %var = 100 [ average are shown. 





a v 


a s 






k v 


k s 


fp 


&c,exc 


Ojfc 




set H 


15.647 


20.610 


1.198 


-1.002 


2.079 


4.451 


1.923 


0.296 


-4.738 


30.717 


seti 


15.834 


21.074 


1.203 


-0.905 


1.305 


-0.812 


1.265 


0.304 


-5.910 


-12.666 


set 2 


15.557 


20.860 


1.213 


-1.037 


1.954 


4.339 


2.583 


0.418 


-4.719 


32.562 


set 3 


15.673 


19.718 


1.234 


-1.005 


1.734 


2.072 


2.071 


0.394 


-1.248 


21.077 


set4 


15.671 


22.448 


1.227 


-0.726 


2.720 


9.059 


2.936 


0.661 


-7.069 


58.552 


sets 


15.577 


21.060 


1.196 


-0.873 


2.398 


6.963 


2.259 


0.334 


-6.359 


44.965 


set6 


15.611 


20.997 


1.158 


-0.963 


2.403 


6.781 


1.457 


0.084 


-8.514 


33.205 


set 7 


15.806 


20.834 


1.204 


-0.962 


1.393 


-0.390 


1.393 


0.289 


-4.988 


-12.233 


set sem j 


15.363 


18.242 


1.221 


-0.904 


1.418 


0.666 


2.577 


0.211 


0.538 


-0.408 


average 


15.638 


20.649 


1.206 


-0.931 


1.934 


3.681 


2.052 


0.332 


-4.779 


21.752 


a 


0.132 


1.077 


0.021 


0.089 


0.479 


3.313 


0.558 


0.149 


2.662 


23.694 


% var 


0.8 


5.2 


1.7 


9.5 


24.8 


90.0 


27.2 


44.9 


55.7 


108.7 



side, along columns 3 and 4, containing the more spherical nuclei, the largest RMS are found 
employing the sets 1 and 7, and in some cases 2 and 6. 

It is tempting to conclude that the LDM is best suited for the description of prolate 
quadrupole deformed nuclei. In order to find support for this conclusion, it is worth to analyze 
the results presented in Tables H and EJ The global RMS for LDM2 is 2.37 MeV, and for LDM3 
is 2.42 MeV. The three fits are pretty close to each other. The smallest RMS are always found 
for the 364 nuclei belonging to group 7, the more prolate deformed, with 0.65 MeV and 0.63 
MeV for LDM2 and LDM3, respectively. The masses of the 332 spherical nuclei included in 
group 3, and the 185 semi-magic nuclei can hardly be fitted with and RMS smaller that 2.0 
MeV. Also the correlations between groups are similar for these two other models. 

This is the most relevant result reported in this contribution: the Liquid Drop Model 
is best suited to describe the masses of prolate deformed nuclei than of spherical 
nuclei. 

3.5. Comparison of the three LDM formulas 

While being close to each other, the three LDM formulas employed in this work have differences 
which are worth to be studied in detail. A statistical analysis of the parameters of each model 
probed to be useful in previous studies |llj . From the nine values of each parameter, their 
average, fluctuations and percentage fluctuations are presented in the last three rows of Tables 
ElElandEl 

The LDM1, Eq. (TjQ) and LDM2, Eq. ([3]) are very similar. Their volume, surface and Coulomb 
parameters, listed in Tables [3] and [5j have nearly equal numerical values, and are all of them very 
stable, with a dispersion along the nine sets smaller than 2%. The pairing parameter fluctuate 
between 10% and 13%, and the coefficient of the symmetry term around 6-7%. The surface 
symmetry term y in LDM1 is the most unstable, with fluctuations of 22%, while its counterpart 
in LDM2, k, has fluctuations of 12%. In this subtle sense, LDM2 could be considered more 
stable than LDM1. It could be useful to perform a deeper comparison of the surface symmetry 
terms, which in one model was obtained asking for consistency between the volume and surface 



contributions, while in the other was designed to incorporate the isospin dependence explicitly. 

The LDM3, Eq. Q, has also stable volume and Coulomb (charge radius) terms, but the 
surface term has fluctuations of the order of 5%, as seen in Table The pairing parameter, 
as in the other two models, fluctuates around 10%. On the other hand, the model has ten 
parameters, and the remaining six have enormous fluctuations in both magnitude and sign. 
These instabilities of the model parameters could be interpreted as a weakness which should be 
addressed. 



4. Shell effects 

The main obstacle for an accurate description of spherical nuclei employing a Liquid Drop Mass 
formula are the shell effects around closed shells. 

In the literature many different ways of implementing shell corrections to the LDM can be 
found; in general, these methods are rather laborious. A simple method was proposed in Refs. 
|15t 116] based on counting the number of valence nucleons. This shell correction is linear and 
quadratic in the total number of valence nucleons n and z, 

Eldmm = E LD Mi + h(n + z) + b 2 {n + z) 2 (8) 

where n and z are the numbers of valence neutrons and protons (particle- or hole-like) and 6, 
are parameters. Inclusion of these two terms in the LDM mass formula (HJ reduces the rms 
deviation from 2.39 to l.OSMeV. 

Following Ref . [13] , we employ also an upgraded version of the terms © , which is suggested 
by the microscopic mass formula of Duflo and Zuker |10t [22] : 

Eldmm 1 = Eldmm + aiSt + a2(S , 2) 2 + a^Ss + a np S np (9) 

where 

nn zz nn(n — n) zz(z — z) nn zz 

2= D~ + D~ J 3 = D + D ' np= D~D~' ( > 

^z - L -'n ^z ^z 

with n = D n — n and z = D z — z, where D n (D z ) is the degeneracy of the neutron (proton) 
valence shell. 

They include 2-, 3- and 4-body terms. The quadratic term is associated to configuration 
mixing and the cubic one to a genuine three body force [18] . 

In Table [8] the RMS of the best fits for all the nuclei, for the seven sets of nuclei grouped 
according to their deformations, and for the semi-magic nuclei, are presented for the Liquid 
Drop Model, Eq. ([JJ, for the Modified Liquid Drop Model, Eq. ([S]) and for the model including 
3- and 4-body terms, Eq. The first row corresponds to the RMS listed in Table [2] in bold 
face numbers. It is clear that the inclusion of microscopic terms improves the fits. The global 
RMS, for all nuclei, diminished form its LDM value of 2.39 MeV to 1.07 MeV and 0.89 MeV. 
The most impressive reductions in the RMS are found in the spherical nuclei grouped in region 
3, which drops from 2.06 MeV to 1.01 MeV and 0.90 MeV, and for the semi-magic nuclei, whose 
RMS diminishes from 2.11 MeV to 1.04 MeV and 0.82 MeV. 

On the other hand, the nuclei in region 7, the most prolate deformed, are the best fitted in the 
LDMM, while the spherical nuclei in region 3 and the semi-magic nuclei have the largest RMS. 
The inclusion of 2-, 3- and 4-body terms in LDMM' seems to succeed in introducing deformation 
effects. Regions 1 and 4 to 7 have all RMS between 562 and 623 keV. Spherical and semi-magic 
nuclei remain to be those with the largest RMS. 



Table 8. RMS (in keV) of the best fit for each of the nine groups (columns), employing Eqs. 
©, © and ([SD, respectively. 





all 


1 


2 


3 


4 


5 


6 


7 


semi-magic 


LDM1 


2387 


1313 


1676 


2063 


1746 


1053 


869 


746 


2113 


LDMM 


1075 


797 


962 


1007 


828 


711 


792 


616 


1038 


LDMM' 


888 


623 


741 


902 


634 


562 


620 


575 


817 



5. Conclusions 

Along this contribution we have shown that the Liquid Drop Model is best suited to describe the 
masses of prolate deformed nuclei than of spherical nuclei. The analysis was performed employing 
three different Liquid Drop Mass formulas. With them, the are nuclear masses nuclei grouped 
in eight sets with similar quadrupole deformations were fitted. For the three LDM models it was 
found that the masses of prolate deformed nuclei can be described with remarkable precision 
for a LDM, with an RMS smaller than 750 keV, while the masses of spherical and semi-magic 
nuclei are those worst described, with RMS larger than 2000 keV. 

The dispersion of the parameters of the three models were studied comparing the fits for the 
different groups of nuclei. We found that in the three the surface symmetry term is the one 
which varies the most from one group of nuclei to another. In the model of Ref. (17] . isospin 
dependent terms were found to exhibit strong changes, making this model the least robust of 
the three under this criterion. 

The inclusion of shell effects allows for better fits, which continue to be better in the prolate 
deformed nuclei region. The Duflo-Zuker model is based in a microscopic description of shell 
effects, and describes deformation through a change in valence occupations. It remains a 
challenge to see if the DZ mechanism to incorporate deformation effects can be successfully 
employed by other models. 
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